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20.  ABSTRACT  (Continued) 

stable  or  unstable.  These  characteristics  are  demonstrated  by  examining  three  general  types  of 
perturbations.  We  further  show  that  deformations  of  the  vortex  sheet  or  the  fluid  layer  affect 
the  flow  field  in  two  ways:  disturbing  the  total  pressure  field  and  perturbing  the  centrifugal 
force  field  created  by  the  azimuthal  components  of  the  velocity  and  magnetic  fields.  The  latter, 
even  though  it  seems  to  be  straightforward,  has  often  been  overlooked  in  previous  analyses. 

It  is  shown  that  failure  to  consider  such  a  perturbation  to  a  stable  centrigufal  force  field  will  lead 
to  the  improper  destabilization  of  modes  with  smaller  axial  and  azimuthal  wave  numbers. 


SaeuAITV  CLASSIFICATION  OF  THIS  FAOS (*•>•»  D—  Mnttod) 


ii 


CONTENTS 


INTRODUCTION  . 

GOVERNING  EQUATIONS  AND  INTERFACIAL  CONDITIONS 

A  GENERAL  TYPE  OF  VORTEX  SHEET  . 

CONCLUSION . 

REFERENCES  . 


ale-  Per 

■  (x:\kl 
'  T  ‘  » 

■nr  'ce<5 


:  r.rifcutlon/ _ 

.  ..tilaMlity  Codes 
[Avail  and/or 
Dist  i  Special 


STABILITY  CHARACTERISTICS  FOR  FLOWS 
OF  THE  VORTEX  SHEET  TYPE 


INTRODUCTION 

The  density  variations  of  stratified  fluids  play  a  subtle  but  crucial  role  in  flow  characteristics.  In 
stratified  shear  flows  under  gravitational  influences,  for  example,  the  density  variations  enter  the  prob¬ 
lem  through  their  interaction  with  the  vertical  displacements  in  the  flows.  The  density  variation,  usu¬ 
ally  very  small  in  atmospheric  and  oceanographic  studies,  is  the  basis  for  stability  conditions  like  the 
Richardson  criterion.  Such  a  density  variation,  no  matter  how  small,  cannot  be  ignored  because  of  its 
subtle  interaction  with  the  perturbation  in  the  vertical  direction. 

In  the  case  of  vortex  motions,  a  centrifugal  force  field  which  depends  on  the  radial  position  will 
be  induced  by  the  rotation  of  fluids.  Because  of  the  complication  of  the  geometry  and  the  dependence 
of  the  centrifugal  force  field  on  both  rotation  and  density,  the  "buoyancy"  effect  and  the  inertia  effect 
are  sometimes  difficult  to  distinguish.  The  rotation  of  fluid  particles  plays  a  dual  role  in  flow  charac¬ 
teristics.  While  the  rotation,  interacting  with  the  density,  generates  a  force  field  to  supply  radial  "buoy¬ 
ancy"  effects  which  can  either  stabilize  or  destabilize  the  flow,  the  velocity  gradient  creates  shear  effects 
which  always  destabilize  the  flow.  In  other  words,  two  types  of  instability  mechanisms,  the  centrifugal 
one  and  the  shear  one,  are  conveyed  by  the  rotation,  and  they  are  not  as  distinct  as  in  the  case  of  two- 
dimensional  parallel  flows. 

In  an  early  paper  on  interfacial  conditions  of  a  cylindrical  vortex  sheet,  Fung  (1980)  showed  that 
perturbations  to  the  flow  disturb  both  the  pressure  field  and  the  centrifugal  force  field  which  is  created 
by  the  fluid  rotation  and  the  azimuthal  magnetic  field.  The  latter  stabilizes  or  destabilizes  the  flow 
depending  on  whether  the  force  field  generated  is  centrifugally  stable  or  unstable.  The  perturbation  to 
the  centrifugal  force  field  is  therefore  essential  to  the  flow  characteristics  especially  when  rapid  changes 
of  flow  quantities  exist  within  a  thin  layer  of  fluid.  Such  a  perturbation  from  the  centrifugal  force  field 
seems  to  be  straightforward  but  is  sometimes  easily  overlooked.  Sufficient  care  should  be  taken  when 
analyses  involving  discontinuous  quantities  are  performed. 
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In  their  analysis  of  inviscid  instability  of  two-dimensional  vortex  type  flows,  Michalke  &  Timme 
(1967)  failed  to  consider  the  effect  of  perturbations  to  the  centrifugal  force  field  generated  by  the  vor¬ 
tex  motion  and  came  to  a  conclusion  all  modes  were  unstable.  In  a  later  paper  by  Leibovich  (1969)  on 
stability  of  inviscid  rotating  coaxial  jets  in  stratified  fluids,  the  same  perturbation  to  the  centrifugal  force 
field  was  again  omitted.  The  author  found  that  the  stabilities  uncovered  by  the  analysis  were 
anomalous  and  therefore  concluded  that  errors  might  result  if  a  thin  but  stable  fluid  layer  was  replaced 
by  a  vortex  sheet. 

When  discontinuities  of  flow  quantities  exist  in  a  cylindrical  interface,  instabilities  are  likely  to 
occur  because  of  any  unbalanced  centrifugal  forces,  and  the  sharp  velocity  gradient  present  at  the  inter¬ 
face.  As  mentioned  earlier,  though  the  velocity  gradient  always  stabilizes  the  flow,  the  centrifugal  force 
induced  by  vortex  motions  can  either  stabilize  or  destabilized  the  flows.  As  for  the  cases  considered  by 
Michalke  &  Timme  (1967)  and  by  Leibovich  (1969),  the  centrifugal  force  field  did  have  stabilizing 
effects  and  should  have  stabilized  modes  with  smaller  wave  numbers.  This  phenomenon  will  be 
demonstrated  by  examining  a  general  class  of  vortex  sheet  type  flows.  The  stabilizing  or  destabilizing 
effect  of  the  centrifugal  force  field  will  be  revealed  by  the  perturbation  of  the  field  at  the  interface.  A 
centrifugally  stable  force  field  created  by  the  rotation  and  the  azimuthal  magnetic  field  may  not  always 
offset  the  shear  instability  of  the  vortex  sheet,  but  certainly  will  stabilize  disturbances  corresponding  to 
longer  wave  lengths. 

GOVERNING  EQUATIONS  AND  INTERFACIAL  CONDITIONS 

The  governing  stability  equations  for  a  general  class  of  incompressible  vortex  flows  with  radius- 
dependent  density,  velocity  and  magnetic  fields  in  a  cylindrical  coordinate  system  are  given  as: 

(N2  -  NA)D*  [■—  -  —  (NCI  —  NASlA)  -d  -  -j-  jfc2  +  ^\q  (1) 
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Here  N  -  kW  +  mil  -  w  is  the  Doppler-shifted  frequency,  NA  -  kWA  +  milA  the  Alfven  frequency, 
o)  the  complex  eigenfrequency,  p0(r)  the  density  of  the  fluid,  il  the  angular  velocity,  W  the  axial  velo¬ 
city,  H  A  the  angular  Alfven  velocity,  WA  the  axial  Alfven  velocity,  u  the  perturbation  velocity  in  the 
radial  direction,  q  the  perturbation  of  the  total  pressure  (including  the  magnetic  pressure),  k  the  axial 
wave  number,  m  the  azimuthal  wave  number,  T  the  surface  tension  at  the  interface  located  at  a  radial 

position  r  —  i?,  5  the  Dirac  Delta  function,  D  —  and  D*  —  D  +  — .  The  Rayleigh-Synge  discrim- 

dr  r 

inant  is  defined  as 

D[p0(r2n)2] 


and  the  Alfven  discriminant  as 


—D(poa2). 


The  corresponding  boundary  conditions  are  the  perturbations  vanish  at  the  inner  and  outer  boundaries. 

The  interfacial  conditions  for  possible  discontinuities  at  the  vortex  sheet  or  cylindrical  fluid  layer, 
obtained  by  integrating  Eqs.  (1)  and  (2)  across  the  interface  or  by  considering  the  total  force  balance  at 
the  interface,  are  as  follows: 


<  — >  -  0 
N 


H*l.  I 


<p0r(il 2  —  n^)>  +  “TjCk2  +  m2  —  1) 


where  k  «  kR  and  <d>>  —  —  ^U?_0)  denotes  a  possible  jump  condition  at  the  interface. 

Readers  are  referred  to  Fung  (1980)  for  the  detailed  derivation  and  assumption  for  Eqs.  (1)  and  (2), 
and  the  discussion  on  the  two  interfacial  conditions  in  Eqs.  (3)  and  (4). 

For  the  convenience  of  mathematical  operations  and  discussion,  we  define 

Fc  -  <p0r(Cl2  -  il  j)>  +  ~zj(k2  +  m2  -  1)  (5) 


to  denote  the  centrifugal  force  and  surface  tension  effects  at  the  interface.  The  interfacial  condition  (4) 
now  written  as 


3 
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can  then  be  viewed  as  the  dynamical  balance  condition  between  the  perturbation  of  the  total  pressure 
(including  the  magnetic  pressure)  and  the  perturbation  of  the  unbalanced  centrifugal  forces  and  surface 
tension  at  the  interface. 


A  GENERAL  TYPE  OF  VORTEX  SHEET 

The  centrifugal  force  enters  the  dynamic  interfacial  condition  through  its  interaction  with  the 
Lagrangian  displacement  in  the  r-direction.  The  interaction  represents  the  influence  on  flow  stability 
due  to  the  perturbation  to  the  centrifugal  force  held.  To  understand  such  an  influence,  we  will  analyze 
a  general  class  of  vortex  sheet  type  flows  subject  to  different  types  of  perturbations.  Analytical  solu¬ 
tions  for  some  particular  flow  profiles  will  be  obtained  to  verify  the  conclusion  of  the  analysis.  Stability 
domains  will  also  be  discussed.  The  errors  that  result  from  omitting  the  perturbation  to  the  centrifugal 
force  field  in  earlier  analyses  [Michalke  &  Timme  (1967);  Leibovich  (1969)]  will  be  discussed  and 
corrected. 

The  general  type  of  vortex  sheet  profile  to  be  considered  has  two  flow  regions  with  their  steady 
state  interface  located  at  r  —  R.  The  flow  properties  in  the  inner  region  are  all  constant,  i.e., 

P0(r)  -  pi 

ft(r)  -  n, 

W(r)  -  W{  for  0  <  r  <  R 

n^(r)  — 

WW-  WAX 


where  the  quantities  with  numerical  indices  are  constant.  The  flow  properties  in  the  outer  region  are 
arbitrary  functions  of  the  radius.  The  solutions  for  the  perturbation  velocity  ut  and  the  perturbation 
pressure  qx  in  the  inner  region  can  be  obtained  from  Eqs.  (1)  and  (2),  and  the  boundary  condition  at 
the  axis  as 


“i 


A  .V, 


2m  (V|fl  |  —  Na  |fl  AX)  + 


kg\lm(kgxr) 

Im(kg\f) 


lm{kgxr) 


(7) 


qx--iA(N t2  -  )g2 p\Im(kgxr) 


(8) 
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where 


r?  -  1  -  4 


NxCl  i  -  NAXClAX 
Nt  -  Nh 


jV,  -  kWx  +  mCl i  -  <u 
Na{  -  kWAX  +  mClAX 

and  Im(kgxr )  is  the  modified  Bessel  function  of  the  first  kind.  The  prime  denotes  the  total  derivative 
with  respect  to  the  argument  of  the  Bessel  function.  Taking  the  derivative  of  Eq.  (7)  with  respect  to  r , 
one  obtains 


2m(Nx(lx  — 


kg\Im(kg\r)  +  +  k2gl\lm(kg\r) 


The  solutions  in  the  inner  region  given  by  Eqs.  (7),  (8)  and  (10)  will  be  used  to  analyze  the  influence 
of  the  centrifugal  force  field  on  flow  stability  subject  to  three  kinds  of  disturbances  at  the  interface:  an 
axisymmetric  perturbation,  an  aximuthal  perturbation,  and  an  arbitrary  perturbation. 

Case  1:  The  axisymmetric  mode  (m  **  0) 

The  solutions  in  the  inner  region  as  described  by  Eqs.  (7)  and  (10)  for  the  axisymmetric  case 
reduce  to 

ux  -  A  Nxkgxl0{kgxr)  (11) 

D*ux  -  A  Nxklg}l0(kgxr)  (12) 


where 


Nx  -  kWx  -  w 


AT.n,-  kwAXdAX 
Nl  -  k2  W\  I 


For  the  convenience  of  mathematical  operations,  we  will  express  the  dynamical  interfacial  condi¬ 
tion  only  in  terms  of  the  perturbation  velocity.  To  do  this,  we  substitute  (1)  into  (4)  for  m  —  0  and 
obtain 
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<P„rHN2  -  *2^2)i>j^)>  +  k2[-^  Fc  -  0. 


The  governing  stability  equation  obtained  by  combining  Eqs.  (1)  and  (2)  for  the  axisymmetric 


mode  reduces  to  the  form 


■M- 


D\Po{N2-  k2Wj)D4-%  -p0*W-  k2fV2)g2  +  m2-<J>  +  *J  -0  (14) 


where 


g2  -  1  -  4 


nsi  -  kwAnA 

N2  -  k 2  W} 


for  R  <  r  <  oo. 


Multiplying 


Eq.  (14)  by  r|-=j. 


where  the  quantities  with  a  bar  are  the  complex  conjugates,  and  integrat¬ 


ing  the  resultant  equation  over  the  outer  region,  we  obtain,  after  applying  the  boundary  condition  at 
infinity,  the  following  integral  equation 


I  F 

-  /~Po*2«t>  -  40 2  -  *A)  rdr  -  0. 


p0(N2-k2W}) 


D*i  +k2g2jf  rdr 


Combining  Eqs.  (11),  (12),  (13),  and  (15),  we  obtain,  by  using  the  complex  conjugate  of  Eq.  (3), 

Ar2K/0(Kgl){plgi/0(Kg1)(A(i2  -  k2Wj i)  -  K/0(«gi)Fc} 

+  p0(N 2  -  k2  fV})l\D*<f»k I2  4-  k2g2\<ftk\2}  -  J  p0k2(< t>  -  402  -  V  A)\4>k\2rdr  -  0  (16) 


where  the  transformation 


Agx  N 


has  been  applied  for  R  <  r  <  «. 
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It  is  very  difficult  to  observe  the  general  characteristics  of  Eq.  (16)  because  the  arguments  of  the 
Modified  Bessel  functions  involve  the  complex  eigenfrequency  &».  To  observe  the  effects  of  the  centri¬ 
fugal  force  held  and  other  flow  quantities  on  the  stability  of  the  flow  and  especially  of  the  interface,  the 
following  special  profiles  are  considered  to  simplify  the  integral  equation  (16).  Let 

WA  —  0  for  0  <  r  <  °° 


and 


n,  -  0 


which  correspond  to  a  vortex  flow  field  with  no  axial  magnetic  flux  anywhere  and  a  core  which  is  not 
rotating.  Equation  (16)  reduces  to  a  quadratic  form 


+  «•* 


0 


(17> 


where 

a*  -  «o  +  Xj  Qkdr 
bk-aaW{  +  f~WQkdr 

ck  -  a„  W 2  -  IkI0(k)]2Fc  +  J  W2Qkdr  -  fp„(<P  -  'V  A)\4>k\2rdr 
a0  —  k2p \kIq (k)  Iq (k )  ^  0 


and 


Qk  -  Po(\D*d>k\2  +  k2\<t>k\2)r  >  0. 


Solving  for  co,  one  finds  that  stability  (corresponding  to  real  values  of  u)  is  guaranteed  when 


+ 


(«*  +  Xj  £?*<*•)  jx,  Po I*  -  *A)\<f>k\2rdr  +  (k/„(ic)]2Fc 


>  0 


(18) 


where 


8*  -  /<?**/  ^2Qkdr  -  (J  WQkdr)\ 

Certain  characteristics  of  the  flow  can  be  observed  from  Eq.  (18).  The  first  pair  of  curly  brackets  con¬ 
tain  information  on  the  axial  velocities  while  the  second  pair  of  curly  brackets  contain  information  on 
the  centrifugal  forces  at  the  interface  and  in  the  outer  region.  Since  both  a0  and  Qk  are  positive 
definite,  the  first  term  in  the  first  pair  of  curly  brackets  represent  the  axial  velocity  difference  at  the 
interface  between  the  inner  and  the  outer  regions,  always  destabilizing  the  flow.  The  second  term  in 
the  first  pair  of  curly  brackets  contains  information  on  the  axial  velocity  in  the  outer  region.  It  can  be 
easily  seen  from  the  Schwarz  inequality  that 

8k  ^  0 


for  all  values  of  W.  Therefore  we  can  conclude  that  the  presence  of  the  axial  velocity  in  the  outer 
region  always  destabilizes  the  flow  except  for  constant  values  of  axial  flows  where  Sk  =  0.  The  tangen¬ 
tial  shears  at  the  interface  and  in  the  outer  region  are  suppressed  since  the  perturbations  are  allowed 
only  in  the  axial  direction.  The  first  term  in  the  second  curly  bracket  is  the  integral  of  the  Rayleigh- 
Synge  and  the  Alfven  discriminants.  For  vortex  flows  subject  to  axisymmetric  disturbances  and  in  the 
absence  of  the  axial  velocity  and  the  axial  magnetic  field,  the  two  discriminants  constitute  the  general 
Michael  condition  saying  that  the  necessary  and  sufficient  condition  for  stability  is 

<D  -  >  0.  (19) 

Equation  (19)  represents  a  state  of  centrifugal  stability  [Fung  (1980)]  and  the  corresponding  integral  in 
the  second  curly  bracket  conveys  information  on  centrifugal  stability  in  the  outer  region.  The  last  term 
in  the  curly  bracket  carries  the  information  on  the  centrifugal  forces  acting  on  both  sides  of  the  inter¬ 
face,  parallel  to  the  first  term  in  the  same  brackets.  As  pointed  out  by  Fung  (1980)  in  his  derivation  of 
the  dynamic  interfacial  condition,  the  jump  condition  arising  from  the  perturbation  of  the  centrifugal 
force  field  is  the  outcome  of  integrating  the  Raleigh-Synge  and  the  Alfven  discriminants  across  the 
interface.  The  last  term  in  the  second  pair  of  the  curly  brackets  in  Eq.  (18)  can  then  be  viewed  as  the 
integral  representation  of  the  generalized  Michael  condition  at  the  interface.  The  sign  of  Fc,  indicating 
whether  or  not  the  resultant  force  at  the  interface  is  centrifugally  stable,  determines  the  stabilizing  or 
destabilizing  effect  on  the  flow.  For  the  present  case 


*1 
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(20) 


Fc-R[p2aj-  (p2ni2-Pln^)]  +  -^-(K2-  i). 


The  surface  tension  always  stabilizes  the  flow  except  for  very  long  axial  wave  lengths  where  k  <  1 .  It 
can  be  shown  (Fung  1980)  that  the  centrifugal  force  balance  condition  in  Eq.  (4)  can  be  obtained  by 
integrating  Eq.  (19)  across  the  interface.  Equation  (4)  can  therefore  be  viewed  as  the  integral  represen¬ 
tation  of  the  generalized  Michael  condition,  representing  a  centrifugally  stable  condition  at  the  inter¬ 
face.  As  predicted  by  the  generalized  Michael  condition,  the  presence  of  the  magnetic  field  in  the  inner 
region  stabilizes  the  flow  while  that  in  the  outer  region  destabilizes  the  flow.  The  rotational  velocity 
immediately  outside  the  sheet,  always  stabilizes  the  flow. 

It  should  be  pointed  out  that  the  flow  profile  being  considered  in  this  axisymmetric  case  can  be 
reduced  to  the  one  examined  by  Leibovich  (1969)  in  his  study  on  hydrodynamic  stability  of  inviscid 
rotating  jets,  if  the  axial  flow  in  the  outer  region  and  all  the  magnetic  forces  in  the  flow  field  are 
deleted.  Eq.  (18)  with  W  **  A  *0  should  have  been  recovered  had  the  correct  interfacial  condition 
described  by  Eq.  (4)  been  used  in  his  paper.  Failure  to  consider  the  centrifugal  force  jump,  which  con¬ 
tributed  considerable  stabilizing  effects  to  the  flow,  led  the  author  to  conclude  that  the  flow  must  be 
unstable  at  least  to  short  waves  and  possibly  to  all  wavelengths.  While  the  vortex  sheet  type  of  flows  is 
susceptable  to  short  wave  perturbations  because  of  the  strong  shear  effect  present  at  the  interface,  the 
centrifugal  force  jump  as  in  the  case  investigated  by  Leibovich  (1969)  will  certainly  stabilize  those  per¬ 
turbations  with  longer  wavelengths.  This  characteristic  is  clearly  shown  in  Eq.  (18)  and  will  be  sup¬ 
ported  by  an  exact  solution  to  be  given  in  the  following. 

Because  of  the  presence  of  the  centrifugal  term  involving  Fc  in  Eq.  (18),  instabilities  for  large 
axial  wave  number  can  not  immediately  be  concluded.  Therefore,  it  is  necessary  to  obtain  solutions  for 
some  specific  flow  profiles  in  the  outer  region  before  the  detailed  stability  phenomena  predicted  by  Eq. 
(18)  can  be  observed.  We  will  examine  the  following  flow  profile 

p0(r)  -  pj 

W{r)  -  W2 

(Hr)  -•  (l2(R/r)2  for  R  <  r  <  oo 

fl^(r)  — 

where  p2,  fF2,  (l2,  and  at®  constants.  The  perturbation  velocity  u2  in  the  outer  region  obtained 
by  solving  Eq.  (14)  is 

u2  -  B(kW2  -  u)kK0(kr). 
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The  stability  boundary  described  by  Eq.  (18)  in  this  case  is 


-  k2(Wx  - 


WJ1  + 


1  k/„(k) 
Pi  4>(k) 


1  kK0(k) 
Pi  K0(k) 


Fc>  0 


(21) 


where  Fc  is  given  in  Eq.  (20).  Since  the  sum  of  the  terms  within  .i.e  square  brackets  is  positive,  the 
stabilizing  effect,  if  any,  will  come  from  the  centrifugal  force  term  Fc.  Even  though  Eq.  (21)  is  likely 
to  be  violated  for  very  large  axial  wave  numbers,  perturbations  corresponding  to  smaller  axial  wave 
numbers  may  certainly  be  stabilized  by  the  centrifugally  stable  forces  at  the  interface.  As  demonstrated 
by  Eq.  (21),  an  erroneous  conclusion  that  no  axisymmetric  modes  are  stable  can  easily  be  reached  if 
the  perturbation  of  the  centrifugal  force  at  the  interface  is  omitted  as  was  done  by  Leibovich  (1969). 
Figure  1  shows  the  stability  domains  for  the  ratio  between  the  centrifugal  force  and  the  axial  velocity 
difference  as  a  function  of  the  axial  wave  number.  Even  though  the  flow  is  unstable  for  very  large  axial 
wave  numbers,  perturbations  corresponding  to  longer  wavelengths  will  certainly  be  stabilized  by  the 
centrifugal  force  jump  at  the  interface. 

Case  2:  The  azimuthal  mode  (k  —  0) 

The  solutions  in  the  inner  region  governed  by  Eqs.  (7)  and  (10)  for  azimuthal  modes  have  the 
forms 


u,  -  A  Ni  rm~l 
DmU\  -  A  AT,  mrm~ 2 


(22) 

(23) 


where 


jV|  -  mfl t  —  w. 


We  will  again  express  the  dynamic  interfacial  condition  only  in  terms  of  the  perturbation  velocity. 
Combining  Eqs.  (1)  and  (4)  for  k  —  0  yields 


<p0r2(N2  -  m2 fl  j) 


D *U 

N 


(24) 


—  <2mrpo(NTl  -  mfl^)>  +  m2Fc 


0. 
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As  for  the  outer  region,  the  governing  equation  obtained  from  Eqs.  (1)  and  (2)  for  azimuthal 
modes  is 


Z>Jp0r2(iV2  -  m2ni)£>*j~:jj 


p/nriJtpofAffl  —  mOj)]  +  p0m2(N2  —  4>  +  4fl2  —  m2£l\  +  ¥,<)  -jj  —  0 


(25) 


for  R  ^  r  <  00. 


.  Jul  . 

by  H*srl,  ir 
(AfJ 


Multiplying  Eq.  (25)  by  rj-£r|,  integrating  the  resultant  equation  over  the  outer  region,  and  applying 
the  boundary  condition  at  infinity,  we  obtain 


rdr 

f  P 

+  f  2mrND(p0Cl)  -  m2[r£l2{DPo)  +  riXp.nJ))  ~  rdr  -  0. 

R  N 

(26) 

Using  the  complex  conjugate  of  Eq.  (3)  and  substituting  Eqs.  (22),  (23),  and  (25)  into  (26),  we  obtain 


B 

Oi 

2 

-  2  b„ 

<0 

K: 

m 

m 

+  cm  -  0 


(27) 


where 


am  -  mp ,  +  fR  Qmdr 

b„  -  (m  -  Dpifli  +  p2fl2  +  fRnQmdr  +  fR  D(p0Sl)\<t>m\2r2dr 

cm  -  (m  -  2)Pl(n,2  -  n^,)  +  2P2(n22  -  n  j2)  -  fc/r 

+  /R"(n2  -  n \)Qmdr  +  fRD\p0((l2  -  n2))\*J2r2dr 
Qm  “  Po(r2\D*<t»m\2  +  m2\<bj2)r  >  0 


(28) 
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and  the  transformation 


.  _  _JL _JL 

A  Rm  N 

has  been  applied  in  the  region  R  <  r  <  <» .  The  flow  will  be  stable  if 

b2  -ancm>  0.  (29) 

Since  Qm  and  therefore  am  are  positive  definite,  we  can  immediately  conclude  that  again  the  sign  of  Fc, 
representing  the  balanced  or  unbalanced  force  condition  at  the  interface,  determines  the  stabilizing  or 
destabilizing  effect  on  the  flow.  The  stabilility  mechanism  conveyed  by  Eq.  (29)  can  not  be  observed 
directly  in  this  case.  To  further  investigate  this  mechanism,  we  utilize  the  transform 

P„  -  p0[r2|Z>*J2  +  («2  -  l)l*J2)r  >  0  (30) 

and  from  Eq.  (28) 

Q„  -  P*  -  PoD(r2\4>m  I2). 

The  coefficients  in  Eq.  (27)  can  therefore  be  rewritten  as 
am  *  "»Pl  +  Sr  Q"dr 

bm-  (m  -  l)p,ft,  +  jR  FlPmdr  (31) 

cm  -  («  -  l)p, (ft,2  -  Oi,)  -  (m2  -  1)  77/?3  +  /(ft2  —  Sl%)Pmdr. 

Substituting  the  coefficients  in  Eq.  (31)  into  Eq.  (29),  we  And,  after  some  mathematical  manipulations, 
that  stability  of  the  flow  is  guaranteed  if 

x,  +  X2  +  x,  >  0  (32) 

where 

x,  -  -  |(m  -  l)p,/ (ft  -  ft  \)2Pmdr  +  8m 
*2  -  [p2  -  Pi  +  /(£p0)r2|u|2<//j[(m  -  l)p, ft  2  +  /ft2?*,*//-] 
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and 


x3  -  mpi  +  J*(?drj[(mj  -  1  )T/R3  +  (m  -  Dpift^i  +  f(ljlQmdr]  >  0 


Bm  -  fPmdrf  Sl2Pmdr  -  j/n^drj2. 

The  roles  played  by  the  flow  quantities  on  the  stability  mechanism  can  be  observed  by  examining  Eq. 
(32).  From  the  Schwarz  inequality  that 

8„  >  0  (33) 

for  all  values  of  11 ,  the  terms  in  xx  in  Eq.  (32)  always  destabilize  the  flow.  Furthermore,  by  comparing 
the  terms  in  xt  with  those  in  the  first  pair  of  curly  brackets  in  Eq.  (18),  we  can  draw  an  analogy 
between  the  two  and  conclude  that  they  both  convey  the  shear  effect  which  destabilizes  the  flow.  The 
first  term  in  X]  is  the  shear  effect  generated  by  the  velocity  difference  at  the  interface.  The  second 
term  in  x(  represents  the  shear  effect  carried  by  the  tangential  velocity  in  the  outer  region  and  is  similar 
to  the  corresponding  term  in  Eq.  (18)  except  that  the  latter  is  induced  by  the  axial  velocity  instead. 
The  term  x2  in  Eq.  (32)  is  the  contribution  to  stability  by  the  density  variation  in  the  centrifugal  force 
field.  Obviously  the  stabilizing  or  destabilizing  effect  depends  on  the  density  difference  at  the  interface 
and  on  the  density  distribution  in  the  outer  region.  Densities  that  increase  with  radius  always  stabilize 
the  flow  and  vice  versa  as  one  would  intuitively  expect.  The  above  discussion  reveals  that  the  angular 
velocity  in  rotating  flows  plays  a  dual  role  in  flow  stabilities:  producing  a  shear  effect  and  inducing  a 
centrifugal  force  field.  The  term  x3  in  (32)  contains  the  information  on  the  surface  tension  and  on  the 
magnetic  field.  The  surface  tension  always  stabilizes  nonaxisymmetric  perturbations  as  is  well-known. 
The  azimuthal  magnetic  fields  in  both  the  inner  and  outer  region  always  stabilize  the  flow  in  spite  of  the 
details  of  the  magnetic  profile.  This  characteristic  is  also  true  for  arbitrary  flows  if  only  the  perturba¬ 
tions  in  the  azimuthal  direction  are  permitted. 

To  further  illustrate  the  stability  characteristics  described  by  Eq.  (32),  we  consider  a  special  flow 
profile 

Po(r)  -  p2 

O(r)  —  0  2  for  R  <  r  <  <».  (34) 

n,(r)-n,2 
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All  the  quantities  with  numerical  indices  are  constants.  The  solution,  obtained  by  solving  Eq.  (25)  for 
the  flow  profile  given  in  Eq.  (34),  is 


(35) 


(36) 


As  previously  predicted,  the  sign  of  Fc  determines  whether  or  not  the  force  condition  at  the  interface 
stabilizes  the  flow.  For  stability  (<o,  —  0)  the  characteristic  equation  from  (36)  requires  that 

—  ( m 2  —  l)p,p2(n,  —  fl2)2  +  (p2  —  Pi)I(/n  —  Dpifli  +  (m  +  l)p2fl2l 

+  /n(pi  +  p2)[(m  —  1  )p | £1  ji  +  (m  +  1)  p2f 1^2  +  (rrt2  —  0  T/R3)  ^  0.  (37) 

As  previously  discussed,  the  first  term  in  the  above  equation  represents  the  shear  effect  at  the  interface, 
always  destablizing  the  flow.  Because  of  the  uniform  rotation  for  r  >  R,  Eq.  (33)  is  identically  equal  to 
zero,  meaning  that  no  shear  effect  exists  in  the  outer  region.  The  second  term  in  Eq.  (37)  is  the  effect 
of  the  density  difference  experienced  in  the  centrifugal  forces  generated  at  the  interface.  Stabilizing 
effects  correspond  to  larger  density  in  the  outer  region.  The  last  term  in  Eq.  (37)  contains  the  informa¬ 
tion  on  the  surface  tension  and  on  the  azimuthal  magnetic  fields  in  the  inner  and  outer  region,  all 
always  stabilizing  azimuthal  perturbations.  Figures  (2a,b,c)  show  the  stability  domains  for  m  «■  2,  5 
and  30  modes  in  the  absence  of  surface  tension.  Both  the  destabilizing  effect  produced  by  the  shear  at 
the  interface  and  the  stabilizing  effect  induced  by  the  azimuthal  magnetic  fields  increase  as  the  wave 
number  becomes  larger.  For  very  large  m,  Eq.  (37)  for  zero  surface  tension  reduces  to 

—  pip2(n,  —  n2)2  +  (p i  +  p2)  (pifi^i  +  p2fl.*2)  ^  o.  (38) 

In  the  absence  of  the  magnetic  force,  the  flow  is  always  unstable  except  for  uniform  rotation  where 

fti  ft2. 

Case  3:  The  arbitrary  mode 

For  simplicity  consider  fl  i  -  fl  A  (  -  0  in  the  inner  region,  and  the  solution  for  the  perturbation 
velocity  reduced  from  Eq.  (7)  is 

u,  -  ANxkl'm(kr)  (39) 
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Fig.  2a.  Stability  domains  for  vortex  flows  subject  to  the  m  m  2  azimuthal  perturbation  < - 

n^t-n^-o; - n^,/n,-o.5,  n^-o-.  — •  n,,-o.  n^/nj-o.s.  - 

n^/n,  -iWft2- os). 
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The  dynamic  interfacial  condition  expressed  in  terms  of  the  perturbation  velocities  is 


<p„r2(N 2  -  N })  D*  {<-  2 mrp0(N{l  -  NACl A)>  +  (k2  +  m2)Fc)  -  0.  (41) 

* 


For  the  outer  region,  when  we  consider  ClA  WA  -  0,  the  governing  equation  for  stability  is 


D  -t~TT  ~  {2mrD  ~T^ 

m 2  +  k2r2  |  N  J  m2  + 


2*  +  p<  ^ I  i  -o.  (42) 


i2  +  k2r2  Y  N 


Multiply  Eq.  (42)  by  /|-=j  and  integrate  the  resultant  equation  over  the  outer  region  to  obtain 

dAjlT  .  r  {2mrD  _£»£*_ 

[  (  m2  +  k2r2  (  Af  jj^  J *  m2  +  k2r2 

»,?  .  4m2n2  t  u  ,  r  P»r2M2  u  j 

+  p  „  Ar  —  <l>  +  — : - —  }  —  rdr  —  1  — ; 7-7  D  *  —  rdr  —  0. 

m2  +  k2r2  N  J  m2  +  k2r 2  \  N 


,n  n 

+  k2r2 


Substitute  Eqs.  (39),  (40),  and  (41)  into  (43)  and  utilize  the  transforms 


Q-».  '»» 


P  -  ri  Ir'IO*!’  +  (*V  +  m‘  -  I)  l+l'lf  »  0 

m 1  +  klr 


where 


^  -  JL 

*  AN 


has  been  applied  in  the  region  R  <  r  <  <».  From  Eqs.  (44)  and  (45),  it  follows  that 

„  p0D(r2\d>\2) 


The  characteristic  equation  obtained  by  using  the  transformation  from  Eq.  (44)  to  (46)  is 


at o2  —  2 bat  +  c  ™  0 


(47) 


where 


a  -  «  +  SR  Q  dr 

b  -  akWx  +  f~  (kWQ  +  m(lP)dr 

+  /"  ( k2W2Q  +  2 kWmSIP  +  m2il2P)dr 


c  »  ak2(  W}  -  Htf,)  -  (k/^Oc)]2 


II 

R 


m2p2£i\ 


m2  +  k 2 


*  + 


2m2Sl2 
m2  +  k2r2 


k2r2 

m2  +  k2r2 


M2  rdr 


and 

a  —  p\Kl'm(n) Im(k)  >  0. 

The  roie  played  by  the  unbalanced  forces  arising  from  the  discontinuities  at  the  interface  can  be 
observed  immediately.  Since  both  a  and  Q  are  positive  definite,  the  sign  of  Fc  therefore  determines 
the  stability  effect  carried  by  the  forces  acting  at  the  interface.  In  the  present  case 

Fc-p2RSl}  +  TU2  +  m2-  I )/R2. 

The  surface  tension  always  stabilizes  nonaxisymmetric  perturbations  as  is  well-known.  The  centrifugal 
force  term  arising  from  the  difference  in  the  angular  velocity  at  the  interface  always  stabilizes  the  flow. 
This  is  also  the  term  omitted  by  Michalke  &  Timme  (1968)  in  their  stability  analysis  of  two- 
dimensional  vortex  type  flows.  Such  an  omission  led  them  to  an  erroneous  conclusion  that  all  azimu¬ 
thal  modes  were  unstable. 

Solving  the  quadratic  equation  (Eq.  (47))  for  to,  we  conclude  that  the  flow  will  be  stable  if 

y\  +  yi  +  >*3  +  y*  >  0  (48) 


where 


y\  -  "  a  J“  k2(W,  -  W)2  Qdr  -  8W  -  Sn 

y2  -  2a  J~  k(W{  -  W)  mSlPdr  +  2  [j*^  kWQdr  JR  mSlPdr  -  J*  Qdr  fR  kWmdPdr 


(K  I  \Kt  oo 

kI'„(k)  lp2  — f~— j  -  P\I„(k)]  +  f  D 

m  +  k1  R 

+  («+/ Qdr)$^ 


2  m2Sl: 


m 2  +  k2r2 


y  4 


i«+/h{  ak2  W2X  +  [k/»]2 


m2  +  k2r 2 


<2P2^2  ,  T  ,  ,  ,  2  ,, 

— - - r  +  - 7  («  +  m  -  1) 

m2  +  K2  Rk 2 


m2  4-  k2r 2 


\<f>\2rdr 


rWdr\fR 


m2(l2Pdr 


and 


8„  -  /  Qdr  f  k2  W2Qdr  -  (j kWQdr^ 

8n  —  J*  Pdr  J  m2fl2Pdr  —  . 

From  the  Schwarz  inequality,  it  follows  that 


8W  >  0 
and 

8n  ^  0 

for  all  values  of  W  and  n.  Several  stability  characteristics  can  be  observed  from  Eq.  (48)  as  follows. 

The  quantity  carries  the  shear  effects  in  both  axial  and  azimuthal  directions,  always  destabiliz¬ 
ing  the  flow.  The  first  term  in  y\  is  the  axial  shear  generated  by  the  velocity  difference  between  the 
inner  and  the  outer  regions.  The  second  term  8W  and  the  third  term  8n  in  y i  are  respectively  the  shear 
effects  produced  by  the  axial  and  azimuthal  velocity  difference  within  the  outer  region.  The 
corresponding  terms  can  be  found  in  the  case  for  the  axisymmetric  mode  and  for  the  azimuthal  mode. 

The  quantity  y2  is  the  shear  effect  interaction  between  the  axial  and  azimuthal  directions.  The 
first  term  in  y2  is  the  interaction  between  the  inner  and  outer  regions,  while  the  second  term  is  the 
interaction  within  the  outer  region.  All  the  terms  in  y2  can  be  positive  or  negative,  depending  on  the 
signs  of  the  velocities  and  wave  numbers,  and  therefore  can  stabilize  or  destabilize  the  flow.  Such 
dependence  implies  whether  or  not  the  axial  or  azimuthal  shear  reinforces  each  other  and  whether  or 
not  the  direction  of  perturbations  strengthens  the  resultant  shear  effect. 

The  first  term  in  y3  is  the  effect  of  density  variation  at  the  interface  and  in  the  outer  region  in  the 
centrifugal  force  field  created  by  the  rotation  of  the  fluid.  Densities  increasing  radially  outwards 
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stabilize  the  flow.  The  second  term  in  >>3  involves  the  integration  of  the  Rayleigh-Synge  discriminant 
over  the  outer  region.  The  condition  for  centrifugally  stable  profile,  i.e.,  <P  ^  0  is  the  precondition  for 
the  sufficiency  condition  of  stability  for  flows  subject  to  perturbations  in  both  the  axial  and  azimuthal 
directions  (Fung  &  Kurzweg  1975).  Positive  values  of  d>  stabilize  the  flow.* 

The  quantity  y4  carries  the  information  on  the  forces  acting  at  the  interface.  The  presence  of  the 
surface  tension  and  of  the  axial  magnetic  field  in  the  inner  region  always  stabilize  the  flow.  The  term 
involving  p2fl22  in  y *  is  the  centrifugal  force  created  by  the  rotation  of  the  outer  region  at  the  interface 
and  is  always  positive.  Because  of  the  nonrotating  core  considered  in  the  inner  region,  any  rotation 
immediately  outside  the  interface  will  have  stabilizing  effects. 

To  further  illustrate  the  stability  characteristics  described  by  Eq.  (48),  we  consider  the  following 
profile 


Po M  “  P2 

fV(r)  -  W2  for  /?</■<  00. 

n  (r)  -  Cl2(R/r)2 

Here  p2.  W 2  and  fl2  are  constant.  The  solution  for  the  perturbation  velocity  obtained  from  Eq.  (42) 
for  the  present  flow  profile  is 

u2-  B(ktV2  +  m(l2-  o>)  k  Km(kr).  (49) 


Equation  (48)  then  has  the  form 
-PlP 1[Jk(W'|-  W2)  -  mftjl2  + 


kI'„(k) 

kK^U) 

Pj  /«(*) 

Pl  KJk) 

k2  U/2  I*  W 
?'k  ./:<«) 


+  pj(l}  +  —7  (k1  +  m1  —  l) 

A 


>  0. 


(50) 


The  first  term  in  Eq.  (50)  is  the  shear  effect  created  by  the  velocity  difference  at  the  interface.  The 
axial  velocity  difference  always  destabilizes  the  flow.  However,  such  destabilization  interacts  with  the 
radial  shear  effect  generated  by  the  rotation  of  the  outer  region.  Whether  or  not  the  interaction  rein¬ 
forces  the  destabilization  depends  on  the  direction  of  the  axial  and  azimuthal  velocities  and  of  the  axial 
and  azimuthal  perturbations.  Since  kI^(k)/  1„(k)  ^  0  and  kK^(k)/ K„(k)  ^  0,  the  second  term  in  Eq. 
(50)  is  always  positive  and  is  contributed  by  the  axial  magnetic  field  in  the  inner  region,  the  rotation  of 
fluid  in  the  outer  region,  and  the  surface  tension  at  the  interface.  All  these  contributions  stabilize  the 
flow. 
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It  is  the  contribution  from  the  unbalanced  centrifugal  force  p 2n22  resulting  from  the  discontinui¬ 
ties  of  the  velocity  and  density  at  the  interface  neglected  by  Michalke  &  Timme  (1967)  and  by  Leibo- 
vich  (1969)  in  their  analyses  of  two  types  of  vortex  sheet  flows.  Such  negligence  led  them  to  conclude 
that  some  of  the  perturbations  which  should  have  been  stabilized  by  the  unbalanced  centrifugal  forces 
at  the  vortex  sheet  were  unstable. 

CONCLUSION 

The  present  analysis  has  demonstrated  several  characteristics  of  vortex  flows  with  variable  densi¬ 
ties.  Unlike  the  velocity  in  the  two-dimensional  shear  flows,  the  rotation  of  vortex  motions  plays  a  dual 
role  in  flow  stability:  producing  shear  effects  which  always  destabilize  the  flow  and  generating  a  centri¬ 
fugal  force  field  which  can  stabilize  or  destabilize  the  flow.  The  stabilization  or  destabilization  depends 
on  whether  or  not  the  force  field  is  centrifugally  stable.  For  flows  of  the  vortex  sheet  type,  the  centri¬ 
fugal  force  arising  from  the  discontinuities  in  the  rotating  velocity  and  the  azimuthal  magnetic  field  at 
the  vortex  sheet  has  significant  influence  on  flow  stability.  The  resultant  direction  of  the  centrifugal 
force  at  the  interface,  dictated  by  the  sign  of  Fc  in  Eq.  (5),  determines  whether  such  force  stabilizes  or 
destabilizes  the  flow.  As  shown  in  Eq.  (6),  the  forces  acting  at  the  vortex  sheet  interact  with  the  per¬ 
turbation  displacement  of  the  deformed  interface.  Smaller  displacements  correspond  to  larger  wave 
numbers.  For  a  given  centrifugally  stable  flow  profile  on  both  sides  of  a  fluid  layer  a  smaller  thickness 
produces  less  centrifugal  stabilizing  effects  and  at  the  same  time  generates  sharper  velocity  gradients, 
both  destabilizing  the  flow.  This  argument  explains  why  a  flow  of  the  vortex  sheet  type  is  most  suscep¬ 
tible  to  instability  corresponding  to  disturbances  with  large  wave  numbers.  For  perturbations  with 
smaller  wave  numbers,  the  magnitude  of  the  displacement  is  greater  and  the  resultant  interaction  will 
have  more  significant  effects  on  flow  characteristics.  Such  effects  will  stabilize  perturbations 
corresponding  to  smaller  wave  numbers  as  shown  in  the  examples  in  the  present  analysis. 

The  deformation  of  the  vortex  sheet  affects  the  flow  in  two  ways:  disturbing  the  total  pressure 
field  and  perturbing  the  centrifugal  force  field  created  by  the  azimuthal  components  of  the  velocity  and 
the  magnetic  flux.  The  latter,  even  though  it  seems  to  be  straightforward,  is  easily  overlooked  as  in  the 
studies  performed  by  Michalke  &  Timme  (1967)  and  by  Leibovich  (1969).  Failure  to  consider  such  a 
perturbation  to  a  stable  centrifugal  force  at  the  vortex  sheet  can  lead  to  the  erroneous  destabilization  of 
certain  modes  corresponding  to  smaller  axial  and  azimuthal  wave  numbers. 
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